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Summary 

A methodology was developed by Lawrence W. 
Rehfield for the structural analysis of composite 
rotor blades. This coupled-beam analysis is rela- 
tively simple to use compared with the alternative 
analysis techniques. The beam analysis was devel- 
oped for thin-wall single-cell rotor structures and in- 
cludes the effects of elastic coupling achieved through 
unbalanced ply orientation. 

This paper demonstrates the effectiveness of the 
new composite-beam analysis method. This is ac- 
complished through comparison of the results of the 
coupled-beam analysis with those of an established 
baseline analysis technique. The baseline analy- 
sis tool is an MSC/NASTRAN finite-element model 
built up from anisotropic shell elements. Defor- 
mations for three linear static load cases are com- 
pared. These loads are centrifugal force at design 
rotor speed, applied torque, and lift for an ideal ro- 
tor in hover. A D-spar designed to twist under axial 
loading is the subject of the analysis. 

At design rotor speed the finite-element analy- 
sis and the coupled-beam analysis indicate, respec- 
tively, 14° and 15° of twist at the spar tip. The 
finite-element model indicates less twist due to rigid 
boundary conditions and wall thickness considera- 
tions. A similar trend is indicated by the applied- 
torque load case. In the applied-lift load case, 
vertical deflections and twist indicated by both 
analysis methods are essentially the same. Results 
indicate the coupled-beam analysis is well within 
engineering accuracy. 

The results presented demonstrate that moderate 
variations in spar twist can be achieved by varying 
the rotor rotational speed. The analysis also pro- 
vides a new and convenient approach for obtaining 
the extensional, torsional, and bending engineering 
stiffnesses. 

Introduction 

A methodology was developed by Rehfield for 
the structural analysis of composite rotor blades 
(ref. 1). The beam force-deformation relationship is 
controlled by cross-section properties that are easily 
defined line integrals around the rotor spar. This 
coupled-beam analysis is relatively simple to use 
compared with alternative analysis techniques. The 
beam analysis was developed for thin-wall single- 
cell rotor structures. It includes the effects of 
elastic coupling achieved through unbalanced ply 
orientation. 

In light of its potential, applications of the 
coupled-beam analysis are needed. One application 
is to design a rotor blade that would change twist as 


a function of rotor speed. A scale model rotor blade 
described in this report was designed to demonstrate 
this coupling. The spar from that preliminary design 
effort was selected as the subject of the analytical 
study presented here. 

The purpose of this paper is to demonstrate the 
effectiveness of the new composite-beam method- 
ology. This is accomplished by analysis of the 
D-spar through comparison of results from use 
of the coupled-beam analysis and an established 
analysis technique. 

Nomenclature 


By,x 

beam flapwise curvature, in 1 

B z ,x 

beam chordwise curvature, in -1 

Cij 

coupled-beam stiffness terms (from 
ref. 1) 

E\\ 

lamina longitudinal modulus, psi 

E 22 

lamina transverse modulus, psi 

EA 

beam uncoupled extensional stiff- 
ness, lb 

EIc 

beam uncoupled chordwise bending 
stiffness, lb-in 2 

Elf 

beam uncoupled flapwise bending 
stiffness, lb-in 2 

G 12 

lamina shear modulus, psi 

GK 

beam uncoupled torsional stiffness, 
lb-in 2 

L 

elastic beam length, in. 

M x 

beam torsional moment, in-lb 

My 

beam flapwise bending moment, 
in-lb 

M z 

beam chordwise bending moment, 
in-lb 

N x 

beam axial force (centrifugal 
force), lb 

Qy 

beam shear, chordwise, lb 

Qz 

beam shear, flapwise, lb 

Sij 

coupled-beam compliance terms 
(from ref. 1) 

U x 

beam extension strain 

y x 

beam chordwise shear strain 

w, x 

beam flapwise shear strain 

w 

distributed beam load, lb/in. 

z tip 

tip deflection due to beam load, in. 


v i2 lamina primary Poisson ratio 

(j) :JC beam twist rate, deg/in. 

Approach 

The baseline analysis tool selected is 
MSC/NASTRAN (refs. 2 to 4). The NASTRAN 
finite-element model is built up from anisotropic shell 
elements. Deformations for three linear static load 
cases are compared. These loads are centrifugal 
force, ideal lift, and applied torque. 

The model rotor blade is shown in figure 1. The 
spar from this design (fig. 2) is the subject of the 
analysis. The finite-element model and the beam- 
element model are developed from the same rotor 
geometry. The spar has a 35.23-in. radius. The 
analytical model is rigid to radial station 5.23 in. 
and has a constant cross section from station 5.23 in. 
to 35.23 in. The model is divided into 50 spanwise 
segments (0.60 in. each). The cross section is defined 
by 16 nodes. 

For comparison purposes the spar is considered 
to be a cantilever and is analyzed as a linear static 
structure. It is essentially rigid from the center of 
rotation to the spar root. In the finite-element model 
a rigid element acts as an end rib connecting the 
outboard end of the beam element to quadrilateral 
elements at station 5.23 in. 

The spar is a high strain graphite and tough- 
ened epoxy composite made of Hercules IM6 fiber 
with Ciba-Geigy R6376 resin. It has a 6-ply layup 
[+20/— 70/+20/-702/+20], with 0° oriented to the 
pitch axis. The cured ply thickness is 0.0055 in. 
The orthotropic material properties are shown in 
table I. 

The beam force-deformation relationship accord- 
ing to reference 1 is 


\N X 1 


U x I 

Qy 


V,X 

Qz 

M x 

~ \Pij\ 

w, x 

4>,x 

My 



M z \ 


- &z,x - 


where \C{j\ is a 6 by 6 beam stiffness matrix. The 
terms in the stiffness matrix are evaluated for the 
D-spar (table II). The nonzero coupling terms (off- 
diagonal terms) are C14, C25, and C36. Thus, cou- 
pling exists between extension and twist, between 
flapwise shear and chordwise bending, and between 
chordwise shear and flapwise bending. Equation (1) 
can then be written as 

N x ~ C\\U, X + C*i4 lb ( 2 ) 


Qy — C22VZ + C25 By,x lb 

( 3 ) 

Qz = C33W x + CsqBz^x lb 

( 4 ) 

M x = CuU, x 4 " in-lb 

( 5 ) 

My = C^V^x + C<y*)By tX in-lb 

( 6 ) 

M z = + C m B z , x in-lb 

( 7 ) 


In order to apply forces and to calculate beam 
deformations, it is necessary to invert equation (1). 
This inversion yields the compliance relationship 


U, x 1 


\N X 1 

K* 


Qy 


= i%] 

Qz 

M x 



My 

- ^z,x - 


M z \ 


where [S^-] — For this particular case, the 

inversion can be accomplished if equations (2) to (7) 
are broken into three pairs of equations with three 
unknowns each. The following expressions for the 
nonzero compliance terms result: 
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Substituting the table II values into equations (9) 
yields the beam compliance values given in table III. 
The beam cross-section force deformation can then 
be written in terms of beam strains and curvatures 
as a function of applied loads: 


U t x — Si\N x + Si4 M x 

(10) 

V x — S-22 Qy + S2*>M y 

(11) 

W,x = S 33 Q Z + S 36 Mz 

(12) 

4 >,X = S 14 Vr + S44MX deg/in. 

( 13 ) 

~ SvsQy d - S bb My in 

(14) 

B z ,x = S 3 eQ z + SqqM z in -1 

(15) 


Application of Analysis 

Centrifugal Load Case 

The model rotor (fig. 1) has a design rotor speed 
of 2077 rpm and weighs 0.0123 lb/in. of span. These 
values are used to calculate the centrifugal force N x 
applied to the spar: 

N x = 934.4 - 0.75287x 2 lb (16) 

where x is measured from station 5.23 in. 

The material density in the finite-element model 
is selected to produce the same weight per unit length 
as the model rotor (0.0123 lb/in.). This is verified 
by comparing the root centrifugal force calculated 
by the finite-element model with the root centrifugal 
force calculated by equation (16). 

Lifting Load Case 

The ideal rotor in hover has a triangular lift 
distribution from radial station 5.23 in. to 35.23 in. 
The total magnitude of the lift load is selected to 
be 10 lb. The shears and moments produced by the 
triangular lift and the blade weight are calculated by 
linear statics. The load w due to lift and blade weight 
can be written in terms of x as 

w = 0.02222x - 0.0123 lb/in. (17) 

Integrating once gives the shear 

Q z = O.Ollllx 2 - 0.0123x - 9.631 lb (18) 
Integrating again gives the moment 

M y = 0. 003704a; 3 - 0.00615x 2 

- 9.631x + 194.5 in-lb (19) 


In the case of the finite-element model, the lift 
is applied as concentrated forces. The forces are 
applied vertically at the spanwise row of grid points 
along the upper surface and at the quarter-chord. 

Torsional Load Case 

The torsional load (torque) is applied at the tip 
of the spar. The torque is constant with respect to 
length and is selected to be 

M x = 50.0 in-lb (20) 

For the finite-element model, the torque is divided 
equally among the nodes at the tip and is applied 
as concentrated moments. An end rib composed of 
shear elements is used in this load case to prevent 
excessive distortion of the cross section. These ele- 
ments provide no resistance to warping in the axial 
direction. 

Calculation of Deflections for Coupled-Beam 
Analysis 

The internal beam forces (centrifugal forces, 
shears, and moments) are calculated for the spanwise 
stations by means of equations (16) to (20) and av- 
eraged over the segment length. The forces are then 
multiplied by the beam compliance matrix to yield 
beam deformation derivatives. These derivatives are 
integrated along the span to produce deflections. 

Results 

Centrifugal Load Case 

In the centrifugal load case, the coupled-beam- 
element analysis and the finite-element model agree 
(fig. 3). The coupled-beam analysis predicts —15.29° 
of twist at the tip, and the finite-element model 
predicts —14.13° at the tip. Figure 4 depicts plotter 
output from the finite-element analysis showing this 
twist. 

The difference in the twist indicated by the two 
analyses can be attributed to several effects. The 
boundary effect induced by the rigid inboard end rib 
has a stiffening effect on the finite-element model. 
The rigid element does not permit any warping of 
the cross section out of its plane. This locally stiffens 
the structure in torsion. The coupled-beam-element 
model used in this case assumes free warping at the 
root. Additionally, the rigid element does not permit 
Poisson contraction at the root. This has a local 
stiffening effect in extension, and thus reduces the 
twist associated with extension. The coupled-beam 
analysis assumes free Poisson contraction. These 
local root stiffening effects are apparent in figure 5. 
The effects begin to die out at approximately radial 
station 10 in. 
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The coupled response of the spar, outboard of 
radial station 10 in., shows good agreement between 
the two analyses. (See fig. 5.) The constant offset 
between the two curves can most likely be attributed 
to wall thickness considerations. The coupled-beam 
analysis is a thin-wall analysis that does not include 
the composite-plate torsional stiffness. The finite- 
element model includes this stiffness. 

Lifting Load Case 

The two analyses also agree well on vertical de- 
flections due to combined lift and blade weight. The 
coupled-beam analysis predicts 6.70 in. of deflection 
at the spar tip. The finite-element analysis predicts 
6.48 in. of deflection at the tip. (See fig. 6.) The 
deflection predicted by the coupled-beam analysis 
is greater than that predicted by the finite-element 
analysis. 

The coupled-beam analysis indicates that spar 
twist is uncoupled from both shear and bending 
(fig. 7). The finite-element analysis indicates a small 
coupling that can be considered as insignificant for 
design purposes (fig. 7). 

Torsional Load Case 

In the applied-torque case agreement is also good. 
Through use of equations (13) and (20), twist <p at 
the spar tip can easily be expressed for N x = 0 as 
follows: 

(j) = S , 4 4 M x L( 180 / 7 r) 

= (0.2003 x 1(T 3 ) (50.0) (30.0) (180 /tt) 

= 17.2° 

The twist predicted by the finite-element model 
is distorted at the tip because of the concentrated 
moments at the end nodes. (See fig. 8.) The apparent 
excellent agreement of the two analyses at the tip is 
due to this distortion. The true agreement between 
the methods is better shown in the twist rate curve. 
(See fig. 9.) Again, the constant offset between the 
two curves in figure 9 can most likely be attributed 
to wall thickness considerations. 

Comparison to Engineering Beam Theory 

Bernoulli beam theory gives us the uncoupled 
relationship between beam curvature bending 
moment M, and bending stiffness El as 

B x = M/EI 

A similar statement can be made for the compos- 
ite beam if coupling effects are ignored. If the 


beam chordwise shear strain is assumed to be zero, 
equation (6) becomes 


Then C55 corresponds to the beam uncoupled flap- 
wise bending stiffness Elf. Alternately, if the chord- 
wise shear force is taken to be zero, equation (14) 
becomes 

By,x ~ 

In this instance, the inverse of S55 corresponds to 
Elf . In the case of the cantilever beam, beam 
strains are not constrained and bending deflections 
due to applied beam loads are sought. For this 
case it is correct to take the beam stiffness from the 
compliance term (S55) rather than from the stiffness 
matrix term (C55). The significance of selecting the 
correct term is illustrated by the particular case. 

For a uniform beam with a triangular load, the 
deflection at the tip (ref. 5) is given as 

Zti P = 11 w tip L 4 /n0EIf in. 

where w tip is the maximum height of the triangular 
lift distribution. For the beam weight alone, the tip 
deflection is (ref. 5) 

z ti p = -wL 4 /SEIf in. 

Superimposing the two cases and considering the 
10-lb triangular lift load and given blade weight of 
0.0123 lb/in. results in 

2f ip = 48 254.6/E Ij in. (21) 

In our particular case, C55 is almost twice as large 
as S^ 1 . Since it is not a realistic constraint to force 
the chordwise shear strain to zero, the beam stiffness 
should not be taken to be C55. For our case the 
chordwise shear force is zero, so Elf is taken to be 

S h5‘ Substituting this value into equation (21) yields 
z tlp = 6.544 in. 

This value agrees well with the deflection predicted 
by Bernoulli beam theory and by finite-element anal- 
ysis. The difference is primarily due to rounding 
off of the constant in the tabulated solution of ref- 
erence 5. Taking the flapwise bending stiffness to 
be C55 would incorrectly indicate that the beam is 
twice as stiff. Similar reasoning shows that since the 
flapwise shear strain is not constrained, the beam 
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chordwise bending stiffness EI C is taken from the 
compliance relationship as S ^ 1 . 

Since beam extension is not constrained, the tor- 
sional stiffness GK is taken from the compliance re- 
lationship as S 44 . Similarly, since beam twist is not 
constrained, the extensional stiffness EA is S^. 

The engineering stiffness constants are summa- 
rized in table IV. It should be emphasized that these 
stiffnesses are, in general, dependent on boundary 
conditions and are therefore not simply cross-section 
properties. 

Conclusions 

Results from the coupled-beam theory agree 
well with the finite-element analysis predictions. 
At design rotor speed the finite-element analysis 
and the coupled-beam analysis indicate, respec- 
tively, 14° and 15° of twist at the spar tip. The 
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Table I. Orthotropic Material Properties 


E\ i , psi 23.1 x 10 6 

£ 2 2,psi 1.4 xlO 6 

2 0.338 

G 12 ,psi 0.73x10° 


Table II. Nonzero Beam Stiffness Terms 


Cl i,lb 0.8332 x 10 ( 


0.1651 x 10 l 


C33, lb 0.3071 x 10 ! 

C 44 , lb-in 2 0.9747 x 10' 

C 55 , lb-in 2 0.1337 x I0 l 

C 66 , lb-in 2 0.1128 x 10‘ 

C14, lb-in -0.6294 x 10 5 

C25, lb-in 0.3147 x 10 5 

C3G, lb-in -0.3147 x 10 r 


Table III. Nonzero Beam Compliance Terms 


Sii.lb- 1 0.2345 x 10- 

5 22 , lb" 1 0.1099x 10“ 

533,1b- 1 0.4561 x 10“ 

5 44 , (lb-in 2 ) _1 0.2003 x 10“ 

5 55 , (lb-in 2 ) -1 0.1356x 10“ 

5 66 , (lb-in 2 ) -1 0.1242x 10“ 

5i4, (lb-in.) -1 0.1513x10“ 

5 25 , (lb-in.) -1 - 0.2585 x 10“ 

5 36 , (lb-in.) -1 0.1273 x 10- 


0.1099 x 10“ 
0.4561 x 10“ 


0.2003 x 10' 


0.1356 x 10“ 


0.1513 x 10“ 


- 0.2585 x 10“ 
. 0.1273 x 10" 





1.14 


2.60 

Figure 1. Model rotor cross section. Dimensions are in inches. 


Sta. 

0.00 


Sta. 



Figure 2. Model rotor spar design. Dimensions are in inches. 








Figure 6. Beam deflection due to lift and blade weight. 
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